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Abstract 

We present a leading order QCD analysis of a sample of neutrino induced charged- 
current events with two muons in the final state originating in the lead-scintillating 
fibre calorimeter of the CHORUS detector. The results are based on a sample of 8910 
neutrino and 430 antineutrino induced opposite-sign dimuon events collected during 
the exposure of the detector to the CERN Wide Band Neutrino Beam between 
1995 and 1998. The analysis yields a value of the charm quark mass of m c = 
(1.26 ± 0.16 ± 0.09) GeV/c 2 and a value of the ratio of the strange to non-strange 
sea in the nucleon of k = 0.33±0.05±0.05, improving the results obtained in similar 
analyses by previous experiments. 

Key words: Charm production, neutrino, dimuon 
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Introduction 



In neutrino-nucleon deep-inelastic scattering, events that present two muons in the final 
state are mainly due to the muonic decay of a charmed hadron produced in a neutrino 
charged- current interaction. The charm quark can be produced both through d — > c and 
s — > c weak currents. The s — > c transition dominates in the antineutrino induced dimuon 
events, whilst in neutrino induced interactions the d — > c and s — > c transition have 
comparable contributions, as the large valence d quark content of the nucleon compensates 
for the Cabibbo suppressed d — *> c transition. Dimuon events in neutrino interactions can 
be used to measure the strange quark content of the nucleon and the value of the charm 
quark mass. Moreover, the charm production mechanism is of great importance for testing 
perturbative QCD predictions. 

Neutrino charm production has been investigated with counter experiments like CDHS [lj, 
CCFR 0, CHARM II |3J, NOMAD ^ and NuTeV [5J, as well as bubble chambers 
(BEBC ©) and nuclear emulsion detectors like E531 [7] and CHORUS [815110] . In this 
paper we present the study of a sample of neutrino induced charged-current dimuon events 
produced in the lead-scintillating fibre calorimeter of the CHORUS detector; the collected 
statistics is the second largest to date. The analysis is performed in the framework of the 
leading order QCD formalism and the result will be expressed in terms of the strange con- 
tent of the nucleus k, the charm quark mass m c , the Peterson fragmentation parameter ep 
and the charmed hadron semileptonic branching ratio B^. 



1 The CHORUS experiment 



The CHORUS experiment has been designed for a short-baseline search of neutrino os- 
cillations in the v^-v? channel at relatively high Am 2 [11]. The detector was exposed to 
the Wide Band Neutrino Beam at CERN from 1995 until 1998. Neutrinos were produced 
at the SPS by a proton beam of 450 GeV; the average energy was 27 GeV with a 
contamination of about 6%. The detector is described in detail in Ref. [12]; it consists of a 
nuclear emulsion target, a fibre tracker, an air gap hadron spectrometer, a lead-scintillating 
fibre calorimeter and a muon spectrometer. 

The 112 ton lead-scintillating fibre calorimeter, described in detail in Ref. [IB], consists of 
three sections (EM, HAD1, HAD2) with decreasing granularity along the beam direction, 
for a total of 5.2 interaction lengths. Owing to its high mass, the calorimeter also pro- 
vides an active target for incoming neutrinos. A large number of neutrino charged-current 
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interactions on lead was collected during the four years of data taking. 

A test beam calibration of the calorimeter performances was carried out during the data 
taking and results are reported in Ref. [13]. Two neural net algorithms were developed 
in order to improve the vertex position resolution and the energy resolution for charged- 
current events originating in the calorimeter. Both algorithms have been tuned for the 
purposes of this analysis on Monte-Carlo neutrino induced events which present more 
than one reconstructed muon in the final state. The vertex finding algorithm achieved a 
longitudinal and transverse position resolution of a(V x ) = 4.0 cm and cr(V y>z ) = 3.0 cm, 
respectively, with an average efficiency of finding the correct calorimeter plane of (82 ± 
1)%. The energy resolution was estimated to be a(E v )/E u = 0.39/^/EZ + 0.14 (with E v 
expressed in GeV). 

A muon spectrometer composed of six magnetized iron toroids instrumented with drift 
chambers and scintillators is located further downstream to identify muons and determine 
their momentum. A resolution varying from Ap/p ~ 0.15 at 20 GeV/c to Ap/p ~ 0.19 at 
70 GeV/c was achieved, as measured with test beam muons. 

The trigger system of the CHORUS experiment is described in Ref. [H] . A dedicated muon 
trigger was set-up to collect a large statistic of "double track" events used both in dimuon 
and trimuon analyses [15]. A total of 4.9 x 10 19 protons were accumulated on the Be target 
in the period 1995-1998 and the overall CHORUS data collection efficiency was about 92% 
with a dead time of 5%. 



2 Theoretical framework 

Neutrino interactions with two muons in the final state mainly occur when the weak 
charged- current on a d or s quark produces a charmed quark and the charmed hadron 
subsequently decays with a probability into a muon and other hadrons. The primary 
muon and the decay muon have opposite electric charge. The large value of the charm quark 
mass gives rise to an energy threshold for the process; this is described in the leading order 
QCD formalism by the so called "slow rescaling" mechanism [T6] . 

2.1 Leading order dimuon cross-section 

The leading order (LO) QCD framework of deep-inelastic neutrino scattering (DIS) has 
been adopted to describe charm quark production. This scheme is reported in Ref. [T7] and 
uses the helicity formalism to describe the neutrino charged-current cross-section. It has 
the advantage of treating in a natural way the different mass scales involved in the heavy 
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quark production process. The neutrino charm production cross-section is written as 



d 2 a 
dxdy 



2G 2 F V -^- (\V cd \ 2 d(x) + \V cs \ 2 s(x) 

IT v 



cosh if) " 



+ 



m 2 sinh if) 
2Q 2 2 



(1) 



in terms of the Bjorken variables x and y and the Fermi constant Gp. V c d and V^ a are the 
CKM matrix elements, d(x) and s(x) being the momentum distributions of the scattered 
quarks. The angle if) is related to E u , and Q ,the incoming neutrino energy, the leading 
muon energy and the W boson four-momentum, respectively: 

cosh if) = = (2) 

lQ* + (E v -E,) 2 



This reduces to when m c vanishes. With m q the mass of the scattered quark and M 
the nucleon mass, the scaling variable \ is defined as 

(Q 2 -m 2 q + m 2 c ) + A(-Q 2 im 2 ql ml) 
X = V W2 , (3) 



where 



1 1 I 1 _M 2 
r\ 2x V 4x 2 Q 



z = - + \l— + 7 v ( 4 ) 



and the function A is given by 

A(a, b, c) = (a 2 + b 2 + c 2 - 2(ab + bc + ca)) 1 . (5) 



In the limit M 2 /Q 2 — > the variable x becomes the usual "slow rescaling" variable |18j . 

When dealing with experimentally measured quantities, the expression ([1]) has to be cor- 
rected for the electromagnetic radiative processes arising in the final state whose effect 
is to lower the primary muon energy and, correspondingly, to raise the hadronic shower 
energy. A correction for such an effect has been applied according to the prescription of 
Bardin [THj. 
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2.2 Parton distributions 



To interpret the experimental data in terms of charm quark production described by Eq. ([I]) 
a parametrization of the quark content in the target nuclei is required. The GRV94LO 
parton distribution functions derived from experimental deep-inelastic scattering data has 
been used (20]. To take into account the non-isoscalarity of the lead-scintillating fibre 
calorimeter, the total valence quark content has been parametrized as 

X g v alO, Q 2 ) = ~7~~ X M va l + X d va \ , (6) 



where, due to strong isospin symmetry, the down quark content of the neutron is taken to 
be the same as the up quark content of the proton. The strange quark content in protons 
and neutrons is assumed to be the same. The non-strange sea-quark contents of the nucleon 
is assumed to be symmetric between up and down quarks. The strange quark content of 
the nucleon is described by the parameter 



K 



J [xs(x, nl) + xs(x, ii 2 )] dx 



J xu(x, ill) + xd(x, nl) 



dx 



(7) 



where /i 2 is an arbitrary reference scale chosen to be 20 GeV 2 . In the following, an SU(3) 
flavour symmetric sea is assumed; this means that the strange and the anti-strange content 
of the nucleon are the same and have the same x dependence. 

An additional parameter a is often introduced in similar analyses that allows the strange 
quark to behave differently from the up and down quarks. This is obtained by multiplying 
the strange quark parton distribution function (pdf) by the factor (1 — x) a . In the following 
the weighted average a. — 2 ± 1 of the values given by the CHARM II and CCFR leading 
order analyses is used. 



2. 3 Charm quark fragmentation and meson decays 



The non-perturbative processes that act to "dress" the bare charm quark produced in the 
neutrino interaction are known as fragmentation or hadronization. The fragmentation is 
usually described as a function of the variable z defined as the fraction of the longitudinal 
momentum taken over by the charm hadron h c : 

z = P L (h c )/Pr x (h c ) , (8) 
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where P™ ax (/i c ) is the maximum hadron longitudinal momentum relative to the W-boson 
direction in the boson-nucleon center of mass reference frame. The Peterson parametriza- 
tion of the charm quark fragmentation is written as 



D(^,e P )cx^ 1 (l-i- T ^) 2 , (9) 

where ep is a free parameter to be determined from the data. In leading order QCD, the 
neutrino cross-section can be factorized as 

d 3 a d 3 a LO 

" (uN -> (ihcX) = , 7, {vN -> /icX)D(z) . (10) 



dxdydz ' dxdydz 

This is no longer true in the next-to- leading order formalism [22j . The transverse mo- 
mentum of charmed hadrons with respect to the W-boson direction is on average small 
((Pt) ~ 0.2 GeV 2 /c 2 ) and is assumed to be distributed as dN/dp\ oc e~ 6p T. The parameter 
b is taken to be 1.1 GeV 2 /c 2 following Ref. 



The fragmentation variable z cannot be determined on an event-by-event basis since the 
direction of the charm quark or charmed hadron cannot be measured directly. The visible 
fragmentation variable is thus defined as 



E 



f i2 



(E„ 2 + E, 



11) 



where E^ 2 and E s ^ are the energy of the muon coming from the charmed hadron decay 
and the hadronic shower final state energy, respectively. The relation between z and z v ; s 
has been studied with the aid of Monte-Carlo simulations. 



3 Data sample and event selection 



The dimuon sample was collected during four years of data taking using a dedicated trigger 
setup based on the energy released in the fibre calorimeter and on the presence of two or 
more well separated tracks in the spectrometer [H]. A total of 6.6 x 10 6 events tagged as 
dimuons by the trigger were reduced to about 5.0 x 10 5 requiring that at least two muons 
were well reconstructed in the muon spectrometer. For each event, the reconstructed muon 
energy E^, the angle of the muon with respect to the neutrino beam direction 9^ evaluated 
at the interaction vertex and the hadron shower energy E^ are used to derive the following 
kinematic variables: 

• E™ = E s h + E^i + E^2, the visible neutrino energy 



S 



• Qlis = 2-BJ xs -E / xi(l — cos# M i) — m 2 , the visible four-momentum transfer squared 

• ^vis = Qlis/2M(E sh + E/jq), the visible Bjorken x 

• ^vis = E^/iEgh + Epz), the visible fragmentation variable, as reported in Eq. [TT1 

The final sample of dimuon events was selected imposing the following selection criteria: 

CI: at least two muons must be well reconstructed, i.e. a minimum of five measured track 
points in the drift chambers placed in the gaps between the iron magnets is required for 
both muons; 

C2: the longitudinal position of the event vertex (V x ) determined by the Neural Net algorithm 
must be inside a fiducial volume given by 290.7 cm < V x < 384.0 cm (this corresponds 
mainly to the HAD1 part of the calorimeter); the transverse vertex positions (V y and 
V z ) are required to satisfy —120 cm < V y>z < 120 cm to ensure hadronic shower lateral 
containment; 

C3: the shower energy must be in the range 5 GeV < E sh < 150 GeV; 

C4: to ensure a good reconstruction quality and control of the acceptances, a cut E^ > 5 

GeV when extrapolated to the vertex position in the calorimeter is applied. This cut 

also reduces the meson decay background; 
C5: the reconstructed neutrino energy must be in the range 10 GeV < E v < 240 GeV to 

ensure a good control of the detection efficiencies; 
C6: a cut Q 2 > 3.0 GeV 2 /c 2 is applied in order to exclude regions in which the Monte-Carlo 

simulation is not reliable, as described in the following section; 
C7: the transverse distance between the two muons extrapolated to the interaction vertex x 

coordinate (d^) is requested to be less than 15 cm to reject background due to hadron 

decays. 

The primary muon is assumed to be the one with the highest pr with respect to the neutrino 
beam direction; this leads to a (95.8 ± 1.0)% efficiency to correctly identify the primary 
muon for neutrino induced interactions and (94.0 ± 1.1)% for antineutrino interactions. In 
cases in which there are more than two reconstructed muons, the two most energetic ones 
are used in the analysis. The average charm induced opposite sign dimuon trigger efficiency 
was evaluated by means of Monte-Carlo simulation on a sample of events surviving the 
selection cuts and found to be (91 ± 6)%. The event statistics is reported in Table [1] where 
the top line indicates the charges of the primary and secondary muons respectively. 

The dimuon sample statistics after the selection is therefore 

N~ + = 10218 ; N— = 1441 ; N + ~ = 975 ; N ++ = 46 . (12) 



Figure [T] shows a scatter plot of the selected events in the (p M i, p^) plane. 
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Table 1 

Statistics for opposite-sign and same-sign dimuon events induced by neutrino interactions. The 
cuts are explained in the text. 



Cut 


-+ 




+- 


++ 


total 


CI 


33251 


30750 


5769 


899 


70669 


C2 


17517 


12387 


2144 


203 


32251 


C3 


16894 


11510 


2023 


186 


30613 


C4 


13948 


10078 


1689 


141 


25856 


C5 


13938 


10051 


1688 


140 


25817 


C6 


12579 


8962 


1270 


110 


22921 


C7 


10218 


1441 


975 


46 


12680 



-150 

> 
a 



a. 

100 



50 



-50 



-100 



-150 



-150 -100 -50 50 100 150 

P„ (GeV/c) 

Fig. 1. Scatter plot in the ,p^ 2 ) plane of the dimuon selected events. The sign of the muon 
momentum represents its electric charge. 

4 Monte-Carlo simulation 



A full simulation of the events induced by neutrinos with the production of a charmed 
quark has been performed. The first step of the simulation consists in the implementation 
of the formalism described in section [2j A specific event generator has been developed by 
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the CHORUS collaboration. It implements JETSET [21] for charm quark fragmentation 
and charmed hadron decays following Eq. ([9]). The incoming neutrino beam is described 
according to a GEANT 3.21 [25] based package that describes the interaction of the SPS 
proton beam with the beryllium target. All relevant components of the beam were part 
of the simulation. The parton distribution functions were implemented with the use of 
the PDFLIB [26] package. The GRV94LO [2D] pdfs were adopted in the production of 
the sample of events used in the fit procedure. The CTEQ3L [27J pdfs have been used to 
obtain an estimate of the systematic error on the final result. The value of the CKM mixing 
matrix elements has been taken from the PDG compilation [28J. The maximum allowed 
incoming neutrino energy was _E™ ax = 300 GeV and the kinematical bounds in the event 
generator were Q 2 min = 2.0 GeV 2 /c 2 and lU 2 in = 4.0 GeV 2 in order to avoid quasi-elastic 
and resonance production processes. 

The generator output is then used as input for the CHORUS event simulation program. 
Based on GEANT 3.21 [25], it contains the complete description of the apparatus and pro- 
vides the response of the various instrumented parts of the detector. The events produced 
in this way undergo all the analysis steps performed for the real data. 

4-1 Background evaluation 

The main source of background to the charm induced dimuon events is given by the muonic 
decay of non-charmed hadrons produced either directly or during the shower development of 
a charged-current event. The background level depends on the probability for a tt^ or a K ± 
to decay before it interacts in the calorimeter or in the first part of the muon spectrometer. 
The relative rates of same-sign and opposite-sign events resulting from this background are 
closely related to the probability of producing positive or negative pions and kaons during 
the hadronic shower development. A standard way to estimate its contribution to the 
opposite-sign dimuon sample relies on the admitted hypothesis that the same-sign dimuon 
events belong to the background. For instance, the background to the charm induced 
opposite-sign dimuon events with a leading \i~ is given by 



where the charm to muon decay process has been disabled in the simulation of the neutrino 
induced charged-current events. A total of about 300 events were obtained after selection 
cuts out of 3 x 10 6 neutrino CC interactions for which a complete MonteCarlo simulation 
of the detector response has been performed. The comparison between same-sign neutrino 
dimuon events in data and MonteCarlo is shown in Figure [21 the fairly good agreement 
between data and MonteCarlo validates the use of relation ([131) . Using this sample of events, 
a value of (Nqq /Nqq ) mc = 0.89 ± 0.10 has been established. Once selection cuts have 




(13) 
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Fig. 2. Comparison between reconstructed quantities for same-sign dimuon events in data (points) 
and MC sample (histogram). 

been applied, the number of events with opposite sign dimuons is thus smaller than that 
of same sign muonj^l; this is mainly due to the fact that the magnet polarity of the muon 
spectrometer was chosen to focus negative charged muons. Consequently, the efficiency to 
reconstruct a fi + is smaller than to reconstruct a The majority of background events is 
in the region of small muon energy where the effect is more pronounced. Similarly the value 
(Nqq /Nqq )mc = 0.11 ± 0.02 is obtained and used to evaluate the non-charm background 
to the opposite-sign dimuon sample with a leading fi + . 



5 Results 



The number of observed dimuon events depends on the charm quark mass (m c ) via the slow 
rescaling mechanism, the amount of strange quark sea (k), the fragmentation parameter 
(ep) and on the branching ratio of charm into muon (B^). From the initial sample of Tabled! 
taking into account selection efficiencies and neutrino-antineutrino cross contamination, a 



The ratios between positive and negative pions and kaons produced in the MonteCarlo sample 
of neutrino induced charged current interactions are ir + /ir~ ~ 1.7 and K + /K~ ~ 1.6 for events 
having > 10 GeV and E n)K > 5 GeV 
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Table 2 

Results of the four parameter fit procedure and parameter correlation matrix given by MINUIT. 

Fit result Parameter Correlation Coefficients 





value 


error 


Global 


m c [GeV/c 2 ] 


K 


ep 




m c [GeV/c 2 ] 


1.258 


0.160 


0.74986 


1.000 


0.522 


-0.732 


0.018 


K 


0.326 


0.048 


0.56786 


0.522 


1.000 


-0.534 


-0.018 


ep 


0.0646 


0.0053 


0.75471 


-0.732 


-0.534 


1.000 


0.031 




0.0959 


0.0038 


0.06935 


0.018 


-0.018 


0.031 


1.000 



total of 8910 ± 180 events with a leading /i and 430 ± 60 events with a leading /i + were 
selected. This represents the second largest sample of neutrino induced dimuons to date. 

5. 1 Maximum Likelihood four-parameter fit 

To extract information on the production mechanism of dimuon events, a maximum like- 
lihood fit of the Monte-Carlo events to the data has been performed. 

A Monte-Carlo reference sample has been used to avoid the painful procedure in which the 
likelihood is evaluated by generating a different sample of events for each set of parameters. 
The charm quark mass was set to m£ = 0.5 GeV/c 2 to avoid a possible bias due to the 
vanishing charm production cross-section below the m T c ei value, while the other parameters 
were set to K ref = 0.45 and ep ef = 0.07. All Monte-Carlo "true" information has been 
retained to be able to weight each event when any of the four parameters is changed 
during the fitting procedure. The reference sample consists of 17439 events that survived 
the selection criteria out of 3.5 x 10 5 fully reconstructed generated events. 

An unbinned likelihood fit of the four unknown parameters has been performed on the 
basis of the distributions of the three experimental observables E™, £ vis and z v ; s . The 
probability density functions of the three observables are constructed by assigning to each 
event in the reference sample a weight equal to the ratio of the cross-section evaluated 
with the current set of parameters to the cross-section obtained with the reference set. The 
absolute normalization of the Monte-Carlo sample was obtained by using the observed 
dimuon production rate with respect to the single muon events and fixing the absolute 
value of the cross-section with the same method as reported in Ref. [29]. 

The fit was performed using MINUIT [30J simultaneously on the four unknown parameters 
and the result of the fit procedure is reported in Table Figures [3] and @] show the com- 
parison between the data and the Monte-Carlo samples produced using the best fit values 
reported in Table [2j The result of the fit is very stable and consistent values are obtained 
by fitting separately each parameter keeping the others fixed. 
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Fig. 3. Comparison between data (points) and Monte-Carlo (histogram) evaluated using the 
parameters found by the fit procedure. The small fraction of events originating from pions and 
kaons decaying into muons is also shown. 

5.2 Systematic uncertainties 

Systematic uncertainties on the result of the fit are mainly due to the theoretical model, 
the background description and the neutrino flux normalization. Theoretical uncertainties 
were taken as the maximum deviation of each fitted value from its central value when 
varying the parameters of the model inside their quoted errors and repeating the fit proce- 
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Fig. 4. Comparison between data (points) and Monte-Carlo (histogram) evaluated using the 
parameters found by the fit procedure. The small fraction of events originating from pions and 
kaons decaying into muon is also shown. 

dure. The contribution of the nucleon pdf has been evaluated repeating the fit procedure 
using the CTEQ3L pdf parametrization. The amount of events induced by pion and kaon 
decays, the background subtraction procedure and the total uncertainty on the overall 
event normalization set the limit on the precision of and of the ratio of the dimuon 
to the charged-current cross-sections a^ + /cr cc . The systematic errors due to the fiducial 
volume and the event selection criteria definition have been evaluated by varying these 
cuts inside the resolutions of the corresponding variables. The major contributions to the 
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Table 3 

Major sources of systematic error uncertainties to the evaluated values of m c , k and ep and B t 





Am c [GeV/c 2 ] 


An 


Ae P 




m c [GeV/c 2 ] 


— 


0.02 


0.004 


0.001 


K 


0.02 


— 


0.005 


0.001 


ep 


0.03 


0.02 


— 


0.001 


Event normalization 


— 


— 


— 


0.003 


Background scale 


0.04 


0.02 


0.004 


0.005 


Fit procedure 


0.01 


0.002 


0.001 


0.0005 


E s h scale (5%) 


0.03 


0.013 


0.001 


0.004 


E sh offset (± 2 GeV) 


0.05 


0.011 


0.001 


0.001 


Radiative corrections 


0.02 


0.005 


0.002 


0.0005 


V cd , V cs , nucleon pdf 


0.02 


0.015 


0.001 


0.002 


Fiducial volume 


0.01 


0.02 


0.001 


0.0012 


Selection cuts 


0.02 


0.003 


0.002 


0.0011 


Total 


0.088 


0.046 


0.0084 


0.0078 



total systematic errors on the four fitted parameters are shown in Table [3j 

Including the systematic uncertainties, the result of this leading order analysis can be 
summarized as follows: 

m c = (1.26± 0.16(stat)± 0.09(syst)) GeV/c 2 
k = 0.33± 0.05(stat)± 0.05(syst) (14) 
ep = 0.065 ± 0.005(stat)± 0.009(syst) 
B^ = 0.096 ± 0.004(stat)± 0.008(syst) 



6 Discussion and conclusions 



The results of this analysis are in agreement with those reported by previous experi- 
ments [Tj2"f3 ,4.5]. as shown in Table HI 

The value of -B^lKdl 2 obtained using the results of this analysis can be combined directly 
with other leading order results [1.2.3J following the prescription given in Ref. [31]. An 
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average value of 

{B^\V cd \% = (0.474 ± 0.027) x 10~ 2 (15) 
is obtained. 

Results of this analysis may also be compared with earlier analyses of events originating in 
the nuclear emulsion target of the CHORUS experiment. Good agreement is found with the 
result = 0.085 ± 0.010 obtained after applying the selection E v > 30 GeV/c [10J. The 
use of that selection for the comparison is justified by the small contribution of events below 
30 GeV in the present analysis, as visible in Fig. [3j The value of the Peterson fragmentation 
parameter ep = 0.059 ± 0.013 reported in Ref. [5] for the same definition of z as used here 
(called ep in Ref. |Sj) is in good agreement with the result found in the present analysis. In 
Ref. [9], the value m c = (1.42 ±0.08) GeV/c 2 was reported, assuming k = 0.38 and a = 1. 
The dependence of m c on k and a is given in Ref. [8]. For a = 2, as used here, and fixing k 
to its fitted value 0.33, the value of m c becomes 1.30 ± 0.08 GeV/c 2 in excellent agreement 
with the value found in the present analysis. 

As an additional result, the large sample of collected dimuon events allowed us to make the 
evaluation of the rate of charm induced opposite-sign dimuon events relative to the charged- 
current events as a function of the incoming neutrino energy. An unfolding procedure has 
been used to take into account all detector effects. The correlation between different energy 
bins has been evaluated using a Monte-Carlo simulation. The correction due to the missing 
energy of the outgoing neutrino in the charm decay has been calculated for each energy 
bin using the Monte-Carlo event sample. The result is shown in Fig. [5] where statistical 
and systematic errors have been added in quadrature. It agrees well with the world average 
dimuon rate shown in Fig. 26 of Ref. [31]. 

In summary, the analysis of the second largest sample of neutrino induced dimuon events 
has been performed using a leading order QCD formalism. The slow rescaling mechanism 
gives a good description of the CHORUS data. Results on the charm quark mass m c , the 
Table 4 

Compilation of latest results on neutrino induced opposite-sign dimuon events. 



Experiment 


N2M 




m c [GeV/c 2 ] 


K 




ep 


This analysis 


8910 


430 


1.26 ±0.18 


0.33 ±0.07 


0.096 ± 0.008 


0.065 ± 0.010 


CDHS p] 


9922 


2123 




0.47 ±0.09 


0.084 ±0.014 


[0.02,0.14] 


CCFR [2] 


4503 


632 


1.3 ±0.2 


0.44 ±0.09 


0.109 ±0.010 




CHARM II [3] 


3100 


700 


1.8 ±0.4 


0.39 ±0.09 


0.091 ±0.010 


0.072 ± 0.017 


NOMAD @] 


2714 


115 


1.3 ±0.4 


0.48 ±0.17 


0.095 ± 0.015 


0.08 ± 0.05 


NuTeV [5] 


2280 


655 


1.3 ±0.2 


0.38 ±0.08 


0.101 ±0.012 
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k 0.01 

a. 



0.008 



0.006 



0.004 



0.002 



w 50 100 150 200 250 300 

Neutrino Energy (GeV) 

Fig. 5. Opposite-sign dimuon production relative cross-section. The solid line represents the Mon- 
te-Carlo prediction obtained using the values of Table [2j 

content of strange quark in the nucleon k, the Peterson fragmentation parameter ep and 
the branching fraction of charmed hadrons into muon, B^, agree with the data obtained 
in similar analyses and improve the overall knowledge on the neutrino induced charm 
production and charm decay mechanisms. 
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